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ABSTRACT 


» 
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M 

The  nonlinear  equations  of  motion  for  an  incompressible  simple  fluid, 
occupying  a fixed  bounded  container,  are  formulated  on  the  basis  of  the 

y 

” finite -linear”  viscoelastic  theory  for  materials  with  fading  memory;  this 
formal  boundary-initial  value  problem  is  then  viewed  as  a nonlinear  abstract 
evolution  equation  on  a certain  Hilbert  space.  It  is  shown  that  a linearized 
version  of  this  evolution  equation  is  associated  with  a linear  dynamical 
system  on  this  Hilbert  space,  and  several  stability  and  asymptotic  behavior 
results  for  this  linearized  problem  are  proved  through  the  use  of  Liapunov 
stability  methods.  On  the  assumption  that  the  original  nonlinear  evolution 
equation  also  is  associated  with  some  dynamical  system  on  the  same  space,  it 
is  shown  that  the  rest  condition  of  the  fluid  is  stable  and  all  motions  are 
bounded.  The  Liapunov  function  employed  for  this  purpose  can  be  interpreted 
as  a mechanical  energy  function  for  the  fluid. 

ft 
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1.  Introduction 


In  this  paper  we  study  a boundary- initial  value  problem  describing  the 
motions  of  an  incompressible  simple  fluid  with  fading  memory,  assuming  a 
"finite-linear"  constitutive  equation  as  formulated  by  Coleman  and  Noll  (1961), 
(1964).  Stability  results  for  linearized  versions  of  this  problem  have  been 
obtained  by  Craik  (1968)  and  Joseph  (1974),  who  employed  spectral  analysis 
for  this  purpose;  more  recently,  Slemrod  (1976),  (1978)  has  performed  a 
stability  analysis  for  one  such  linearized  problem  by  using  the  ideas  of 
dynamical  systems  theory.  We  refer  the  reader  to  Slemrod  (1976)  for  a critique 
of  the  spectral  analysis  approach. 

Here  we  are  primarily  interested  in  the  highly  nonlinear  equations  of 
motion  which  result  from  a careful  formulation  of  the  general  problem  described 
above,  employing  only  physically  reasonable  assumptions . Oui  must  important 
assumption  is  that  these  nonlinear  equations  do  lead  to  a dynamical  system 
on  an  appropriate  state  space,  and  we  motivate  this  assumption  by  also  studying 
a linearized  version  of  the  problem.  Although  our  linearized  problem  is 
closely  related  to  the  linear  problem  considered  by  Slemrod  (1976),  (1978),  we 
are  able  to  show  that  certain  of  the  stability  results  for  our  linearized  problem 
do  carry  over  to  the  original  nonlinear  problem,  whereas  those  of  Slemrod  (1976)  , 
(1978),  Craik  (1968),  and  Joseph  (1974)  apparently  do  not.  Our  approach 
here  is  based  entirely  on  the  ideas  of  dynamical  systems  theory. 

In  §2  we  formulate  the  basic  equations  of  motion  for  a simple  incompressible 
fluid,  based  on  the  "finite-linear"  viscoelastic  constitutive  equation  of 
Coleman  and  Noll  (1961).  Assuming  that  for  all  time  t i 0 (but  not  t < 0)  the 
fluid  is  incompressible  and  fills  a fixed  bounded  container,  we  obtain  a 


i 


2 


formal  boundary-initial  value  problem  with  history  dependence  of  possibly 
infinite  duration.  Choosing  a state  space  equipped  with  a "fading  memory" 
norm  (see  Coleman  and  Mizel  (1966)),  we  view  the  formal  problem  as  an  abstract 
evolution  equation  on  this  space;  tractability  of  this  problem  seems  to  be  highly 
dependent  on  selection  of  the  appropriate  state  space. 

As  we  are  unable  to  prove  that  our  nonlinear  evolution  equation  leads 
to  a dynamical  system,  we  consider  a linearized  version  in  §3.  There  we  show 
that  this  linear  evolution  equation  does  lead  to  a dynamical  system  (on  the 
same  space),  and  we  obtain  certain  results  on  stability,  asymptotic  stability, 
and  exponential  decay  of  motions.  Although  our  linear  problem  is  closely 
related  to  those  of  Craik  (1968),  Joseph  (1974),  and  Slemrod  (1976),  (1978), 
we  make  fewer  assumptions  regarding  the  deformation  history.  Finally,  in  §4 
we  assume  that,  in  a certain  sense,  the  original  nonlinear  problem  is  related 
to  a dynamical  system;  we  then  are  able  to  prove  that  the  rest  condition  is 
stable  and  every  motion  is  bounded  in  terms  of  the  initial  state  of  the  fluid. 

Our  results  are  obtained  through  the  use  of  a Liapunov  function  (see  Hale  (1969)) 
which  we  interpret  as  a mechanical  energy  function  for  the  fluid. 
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2.  Formulation  of  the  Problem 


We  consider  an  incompressible  simple  fluid  occupying  a bounded  domain 
H6R3,  with  C ^-smooth  boundary  F,  for  all  time  t 6 R+  = [o,®).  Following  the 
development  of  Noll  (1958)  and  Coleman  and  Noll  (1961),  (1964),  we  formulate 
in  this  section  a set  of  nonlinear  equations  of  motion  based  on  the  "finite- 
linear"  viscoelastic  theory. 

For  this  purpose,  consider  an  arbitrary  fluid  particle  that  has  position 

T)  = (Tl1>Tl2»1V  € Q aC  time  teR+>  letting  X(T;Tl,t)  = (X1,X2,X3)  €R3  and 

p(T;T\,t)  6 R denote  its  position  and  mass  density,  respectively,  at  time  t€R 

(note  that  t ^ 0 but  T may  be  negative) . The  relative  de formation  gradient 

F(t ;”n,t)  is  the  second  order  tensor  whose  components  are  given  by 

f . . CT ; t ) = b . X.  (T;  7],  t) , where  3.  = . It  is  known  that 

J 1 J o M. 


(2.1)  P Ct  ;TJ,  t)  = p (T ; T|,  t)  det  F(T  ; T),  t) 

for  all  (T,T],t)  €R  X Q X ft+. 

Following  Coleman  and  Noll  (1961),  we  denote  by  C(T;1),t)  the  relative 
right  Cauchy-Green  tensor  with  components  c„  = • Here  and  in 

the  sequel  we  employ  the  convention  of  summation  on  repeated  indices.  Under 
the  assumption  of  isotropy  and  homogeneity  (see  Coleman  and  Noll  (1961),  (1964)), 
the  basic  constitutive  equations  of  the  "finite-linear"  theory  of  viscoelasticity 
state  that,  for  a compressible  simple  fluid,  the  components  of  the  stress 
tensor  S(T],  t)  are  given  by 


s^C^.t)  = - -[p(P(t;T),t))  + j h(s,P(t;Tl,t))[ckk(t  - s;T),t)  - 6kk^ds)  6 


ij 


- [ m(s,p(t;'Tl,t))  [c  (t  - s ;TJ,  t)  - 6 ]ds  , (Tl,t)€nxR+  , 


'ij 


ij' 
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where  6.,  denotes  the  Kronecker  delta  and  p,  h,  m are  scalar-valued  material 
lj 

functions  (see  equation  5.18  in  Coleman  and  Noll  (1961)).  Here  we  assume 
the  fluid  to  be  incompressible  for  T£ft  ; hence,  the  density  is  constant  and 
we  set  p(T;T),t)  = 1 for  t€R+,  T|  € Q,  t£R+,  without  loss  of  generality.  Under 
this  assumption  the  above  expression  for  stress  must  be  replaced  by 

'00 

(2.2)  s (T),t)  = - p(T),t)6  - j m(s)[c  (t  - sjTl.t)  - 6 ]ds  , (Tl,t)€QXR+  , 


where  p(T|,t)  is  a constitutively  indeterminate  quantity  and  m : ft  — ft  is  a 
material  function.  We  assume  that  this  function  satisfies  the  following 
conditions : 

a)  m€C1(R+)  fl  £ (R+)  , 

(2.3)  b)  m is  nonnegative  and  nonincreasing  on  R+  , 

c)  m(s)  > 0 for  s€[0,r);  m(s)  = 0 for  s£[0,r)  , 


where  r > 0 may  be  infinite.  This  mild  assumption  follows  from  the  concept 
of  "fading  memory"  as  stated  by  Coleman  and  Noll  (1961). 

We  remark  that  the  validity  of  equation  (2.2)  depends  on  the  assumption 
of  incompressibility  for  t 6ft+,  but  does  not  presume  incompressibility  for 
T < 0.  Furthermore,  we  note  that,  by  equation  (2.1),  the  assumption 


P (T i "H* t ) * 1 for  t €ft+  implies  that 


det  F (T ;T\, t)  - l/p(T;^,t)  , T €ft 


(2.4) 


det  F(T;T],t)  - 1 , , 


for  all  (11, t)  (tlx  R+. 
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We  assume  the  fluid  container  to  be  fixed  in  some  inertial  system,  and  that 
all  particle  positions  are  measured  relative  to  this  inertial  system.  We  also 
assume  the  body  force  (per  unit  mass)  at  ( T|,  t ) to  be  the  gradient  3^q(Tl,t)  of 
a known  external  potential  q : 0 X ft+  -*  R.  Consequently,  the  balance  law 
for  linear  momentum  takes  the  form 


(2.5)  — ~ X.  (T ;T|, t) 

St  1 


T=t 


d.q(Tl.t)  + djS_(T],t) 


(T],t)  6Q  x R+ 


It  is  convenient  to  introduce  some  new  notation;  let 


(2.6) 


v(T],t)  = X(T;Tl,t)  , 

T=  t 

u(s,Tl,t)  = X(t-s;'Tl,t)  - X(t  ;T),  t)  = X(t-s;T],t)  - T)  , 


for  all  (s,T],t)  X QX  ft+,  For  a particle  with  position  T)  6 at  time  t€R+, 

v(Tj,t)  represents  its  velocity  at  time  t and  u(s,T),t)  + T;  its  position  at 
time  t - s S t.  Some  straightforward  manipulations  show  that 


(2.7) 


£ 

3t‘ 


2 X (T;Tl,t)  - v 01,t)d  v (H,t)  + §£  v (H,t>  , 

T=t  ^ -* 


(2.8)  ui(s,11,t)  = - vJ(Tl,t)[6ij  + 3 u^s.H.t)]  - u^s.H.t)  , 


(2.9) 


ui(0,Tl,t)  = 0 , ^ ui(s,Tl,t) 


s=0 


= - v.(11,t)  , 


(2.10)  cij(t  - s ;Tl, t)  - &iJ  = fki(t  - s;H,t)fkj(t  - s;11,t)  - 6^ 

= 3iuk(s,'ll,t)*djUk(s,'fl,t)  + a.u^s.H.t)  + 3jui(s,H,t) 

for  all  (s.Tl.t)  X flX  R+.  Since  the  first  of  (2.9)  implies  3 u^(0,T|,t)  = 0, 


relation  (2.4)  leads  to  the  condition 


6 


det  1 3 .u . (s  , T\,  t)  + 6 . . j 
J i ij 


det  | 3.u.  (0,T1  t)  + 

J l ij 

det I 6 .,  | = 1 , 

ij 


for  all  t^s^O,  T)  E C2;  hence,  we  also  have 

(2.11)  0 = (fl  det|3jUi(s,Tl>  t)  + &_[  J 


s=0 


aJ  Vi(S)1U) 


s=0 


- - a^CTl.t)  , (Tl,t)€QXR+  , 


where  we  have  used  the  second  of  (2.9). 

Upon  collecting  (2.2),  (2.3),  (2.5) -(2. 11) , and  making  a final  physical 
assumption  that  the  velocity  of  the  fluid  is  zero  at  the  boundary  T,  the 
fluid  motion  is  seen  to  be  described  by  the  formal  evolution  equations 


(2.12)  ^ v.(Tl,t)  = - Vj  (7],  t)  ^jvi  ("H,  t)  - 9.p(Tl,t)  + &iq(Tl,t) 


r 

- 3^  j m(s)[aiuk(s,T\,t)  ■a^uk(s,Tl,t)  + ^u^s.H.t) 


+ ajui(s,11,t)]ds  , T) € n , 


(2.13) 


^ u.(s,H,t)  = - v^.(Tl,t)[6_  + a^u.  (s  , Tl,  t)  3 


subject  to 


- ui(s,Tl,t)  , (s,Tl)  €R+  x 0 , 


aiViCT1, t)  “ 0 » T'6n  > 


(2.14) 


ui(0,T1,t)  = 0 - v1(H,t)  + ui(s,Tl,t) 


, Ti  € n , 


s«0 


/± cn, t)  ■*  o , Tier  , 


for  all  t€W  , with  prescribed  initial  data  (v^(T\,0)  .u^s.Tl.O))  satisfying  (2.14), 
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The  unknown  pressure  gradient  in  (2.12)  creates  certain  difficulties,  and 

we  now  employ  a well  known  device  for  removing  this  term  (see  Fujita  and  Kato 

3 

(1964)).  Consider  the  Hilbert  space  (£^(0))  , equipped  with  the  usual  inner 

3 3 

product,  and  let  II  : (Jb^(O))  ->  (X^(0))  be  the  orthogonal  projection  whose 

range  fi(ri)  is  the  closure  in  (^(ft))^  °f  {v€(Co(fi))  |d^v,,(Tl)  = 0,  T]€Q}. 

3 

The  range  of  II  is  orthogonal  to  the  closure  of  the  set  of  elements  of  (£^(Q)) 
which  are  gradients  of  scalar-valued  functions.  Writing  n as  a symmetric 
second  order  tensor  n and  applying  it  to  (2.12),  we  find  that 

(2.15)  v.(Tl.t)  = - nuv.(Tl,t)b.vx(TU) 

Pr 

- n.^j  j m(s)[SAuk(s,Tl,t)*SjUk(s,Tl,t) + 3yix(s,Tl,t)]ds  , Tj 


Once  a solution  (v,u)  of  equations  (2. 13) - (2. 15)  is  found,  the  corresponding 
pressure  gradient  d^p ( Tj,  t ) can  be  recovered  through  equation  (2.12),  or 
through  the  equation  obtained  by  applying  the  projection  I - II  to  (2.12). 

We  wish  to  emphasize  two  points  about  our  physical  assumptions. 

It  should  be  noted  that  conditions  (2.14)  do  not  require  the  fluid  to 
have  always  occupied  the  domain  ft  in  the  past;  i.e.,  for  s > t s 0,  we 
have  not  assumed  that 


(2.16) 


u(s,T],t)  + Tj  € Cl  , T|  € Cl  , 
u(s,T],t)  = o , Tier  . 


Secondly,  conditions  (2.14)  do  not  require  the  fluid  to  have  always  been 
incompressible  in  the  past;  i.e.,  for  s > t ^ o,  we  have  not  assumed  that 


1 


6 

(2.17)  det|a.u.  (s,T),t)  + 6.  ! = 1 , T)€P  . 

It  is  apparent  that  at  least  the  first  of  these  non-assumptions  is  physically 
important . 

We  now  wish  to  put  (2. 13) - (2. 15)  in  the  form  of  an  abstract  evolution 
equation  on  an  appropriate  state  space.  For  this  purpose,  first  consider 

00  3 o 

the  linear  space  S of  pairs  (v,u)  in  (C  (0) ) x U (P  ) where 

° 0 < 3 < r P 


P = {w€C°°([0,r)  X 0)  IwCs.Tl)  = 0 for  (s,Ti)  £[8,r)  x n)  , 


such  that  v€ft(II)  and  u(0,TD  = 0 = v(Tl)  + ^ u(s,Tp 


s=0 


for  all  T|  6Q.  Recalling 


conditions  (2.3)  on  m(s) , where  r > 0 might  be  infinite,  we  define 


(v,u)||3  = [ 1 v (7]) v (TO  + f m(s)d  u (s,Tl)  u.  (s,Tl)ds|  d£I 
1 'V  1 1 k i k i J 


and  we  let  X denote  the  Hilbert  space  obtained  by  ||  • ||^  - completion  of  S. 
Also,  for  x = (v,u)  6 S we  define 

Nx  s (w,z)  , 


(2.18)  w.(Tl)  = - "^(TO^^OD 

• Vj  J!  m(s)[a<jUk(s,T0 -3 juk(s,l\)  + 3 .u^Cs.TDlds  , 
j o J 

zi(s»TD  s (TO  + 3jui(s,T0  j - u.(s,Tl)  , 


and  we  consider  a metric  on  S defined  by  d^(x,x)  = ||x  - x|j^  + ||Nx  - Nx||^,  x,x€S. 
We  denote  by  i(N)  the  d^  - completion  of  S,  and  we  let  N : (fi(N)  CX)  -*  X be  the 
operator  defined  on  £(N)  by  (2.18).  Clearly,  £(N)  is  dense  in  X,  and  the 
completeness  of  (^(N),d^)  implies  that  N is  closed  as  a (nonlinear)  mapping 
from  X to  X. 


4 


Within  this  setting  we  can  replace  (2. 13) - (2. 15)  by  an  abstract  evolution 


equation  on  X,  given  by 


(2.19) 


x(t)  = Nx(t)  a.e.  t £ R+ 

x(0)  = x €X 
o 


This  highly  nonlinear  evolution  equation  is  very  difficult  to  analyze,  even 
as  to  existence  and  uniqueness  of  solutions.  However,  for  x^  = 0,  we  notice 
that  equation  (2.19)  does  admit  as  a solution  the  rest  condition,  x(t)  = 0 
for  all  t€R+.  In  the  following  section  we  study  a linearized  version  of  this 
equation  and  show  that  it  leads  to  a stable  linear  dynamical  system  on  X. 


i'J 


3.  A Linearized  Problem 

We  are  unable  to  prove  that  the  nonlinear  abstract  evolution  equation 
(2.19)  is  related  to  any  dynamical  system  (see  Hale  (1969));  in  §4  we  will  simply 
assume  that  such  a relationship  does  exist.  To  suggest  the  plausibility  of 
this  forthcoming  assumption,  we  now  perform  an  appropriate  linearization 
of  (2,19)  (on  the  same  Hilbert  space  X)  and  show  that  the  resulting  abstract 
evolution  equation  does  lead  to  a (linear)  dynamical  system  on  X.  We  remark 
that  this  Hilbert  space  is  different  from  the  one  used  by  Slemrod  (1976),  (1978) 
for  a similar  but  not  identical  linear  problem. 

Our  first  step  is  to  return  to  the  formal  equations  (2. 13) - (2. 15)  and 
delete  all  nonlinear  terms;  hence,  conditions  (2.14)  are  retained,  while 
(2.15)  and  (2.13)  are  replaced  by 

a rr 

(3.1)  Ft  vi(T^,t:)  = " ^ii3!  J n>(s)aiu£(s’Tl’t)ds  » > 

J 0 3 


(3.2)  u.(s,T),t)  = 


v.  (Tj,  t)  - jr-  u . (s , Tj,  t) 

X OS  X 


(s,Tl)  <=R+  x fi 


for  all  t€ft  . We  continue  to  assume  that  m(s)  satisfies  conditions  (2.3). 

Recalling  the  linear  space  S and  the  Hilbert  space  X defined  in  Section  II, 
we  define,  for  x = (v,u)  €S, 


Ax  = (w,z)  , 

wi(T[)  s - j m(s ) a^u^ (s , TQ ds  , 


z.(s,Tl)  = - v (ID  - -r-  u (s,Tl)  , 

1 1 OS  1 


(3.3) 
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and  we  consider  a norm  on  S defined  by  ||x'||  = ||xj|-^  + ||Ax||^,  x€S.  We  denote 

by  £(A)  the  )1  *]1  - completion  of  S,  and  we  let  A : (£(A)  ^X)  ->  1 be  the  linear 

A 

operator  defined  on  JO  (A)  by  (3.3).  We  see  that  £(A)  is  dense  in  X and  A is 
closed  as  a linear  mapping  from  X to  X. 

The  formal  linear  problem  (3.1),  (3.2),  (2.14)  leads  us  to  consider 
a linear  abstract  evolution  equation  defined  in  the  Hilbert  space  X by 


(3.4) 


x(t)  = Ax(t)  , t € ft+  , 
x(0)  = xo€£(A)cX 


We  wish  to  show  that  (3.4)  is  related  to  a linear  dynamical  system  on  X,  and 

that  the  motions  of  this  dynamical  system  are  (unique)  solutions  of  (3.4) 

for  x € & (A) . 
o 

We  recall  (see  Hale  (1969),  Yosida  (1978))  that  a dynamical  system 

[T(t)}  , on  a metric  space  X,  is  a family  of  continuous  operators 

T(t)  : X -•  X such  that  T(-)x  : ft+  — X is  continuous,  T(0)x  = x,  and 

T(t  + h)x  = T(t)T(h)x,  for  all  t,h  € »i+>  x € X.  The  mapping  T( • ) x : R+  — X and 

the  set  Y(x)  = U T(t)x  are  called,  respectively,  the  motion  and  the  positive 
t 2 0 

orbit  corresponding  to  the  initial  state  x£X.  A subset  CJC-X  is  said  to  be 
positive  invariant  under  IT  (t)  } t £ ^ if  Y(x)  for  every  x € Q.  A motion  T (• ) x 
is  said  to  be  stable  if,  given  any  c > 0,  there  exists  a neighborhood  (x) 
of  radius  6 > 0 such  that  y€tV(x)  implies  T (t)  y € ft  (T  (t)  x)  for  all  t€R+; 

T(*)x  is  asymptotical lv  stable  if  it  is  stable  and  T(t)y  — T(t)x  as  t — 00 
for  all  y in  some  neighborhood  of  x.  Furthermore,  we  recall  that  [t(t))t2^ 
is  called  a linear  dynamical  system  if  X is  a Banach  space  and  T(t)  : X -•  X is 
a linear  operator  for  every  t€ft+;  the  infinitesimal  generator  B : (£(B)CX)  — X 


of  such  a linear  dynamical  system  is  defined  by 
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Bx  = lim  ^ [T(t)x-x]  , 

t\M)  c 

where  £(B)  is  the  set  of  x € X such  that  this  limit  exists.  If  B is  the 

infinitesimal  generator  of  the  linear  dynamical  system  {.T(t)}t2Q  defined 

on  X,  then  for  every  x^G^CB)  the  motion  T(»)xq  is  the  unique  strong  solution 

of  the  equation  x(t)  = Bx(t) , t ^ 0,  for  the  initial  state  x(0)  = x . 

o 

With  this  terminology,  we  can  prove  the  following  result  for  our  linear 
abstract  evolution  equation  (3.4). 

Theorem  3.1:  For  m(s)  satisfying  (2.3),  the  linear  operator  A:  (i(A) CX)  — X 
is  the  infinitesimal  generator  of  a linear  dynamical  system  lT(t)Jt2Q  on  X, 
with  ||T(t)x||^  £ i| xll^  for  all  t€ft+,  x € X. 

Proof:  As  £(A)  is  dense  in  X,  for  the  theorem  to  be  proven  it  is  sufficient 

(see  Yosida  (1978))  to  show  that  -A  is  ||  • |]^  - accretive  and  that  the  range 
ft (I -A)  = X.  Using  the  natural  inner  product  for  X,  and  considering  arbitrary 
x = (u,v)  €£(A),  we  see  that 

(x,Ax)  = - [ v,  (Tl)'TT  .&  f m(s)3  .u.  (s,Tl)ds  dO 

1 iX  J J0  J x 

' n>(s)akui(s»Tl)  •\^vi(Tl)  + §7  ui(s,Tl)3ds  dn 

= - 1*  V.(T1)*S.  [ m(s)3  u.(s,T))ds  dft 

cn  x J Jo  J 

r rr 

' m(s)S.  u (s,Tl)ds  dfl 

on  k l K i 

* 2 h [6kUi(s>11),Vi(s>11)]dsdn 

" 2 InJ0m'<s>*\ui<s»^,dkui(s»1[1)ds  dfi  S°  5 

hence,  -A  is  accretive. 
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As  -A  is  accretive,  it  follows  that  I -A  possesses  a bounded  inverse 
defined  on  fl(I-A).  As  A is  a closed  operator,  it  follows  that  (I  - A)  ^ is 
closed  and  bounded;  hence,  R(I  - A)  = X if  R(I  - A)  is  dense  in  X.  With  this 
in  mind,  let  (w,z)  be  a fixed  but  arbitrary  element  of  the  dense  set  Sc;X, 
and  consider  the  equations 


(3.5) 


v.(TD+",#3.  ^m(s).S.u  (s,T|)ds  = w (Tl)  , 

1 lXjo  jX  1 


U.  (s.Tl)  + v (T\)  + r—  u.  (s , T\)  = z (S.TD 
1 1 os  1 1 


If  these  equations  can  be  shown  to  have  a solution  (v,u)  belonging  to  £(A) 
we  will  have  shown  that  R(I  - A) 3S,  and  by  the  denseness  of  S in  X it  will 
follow  that  ft(I-A)  = X.  Formally,  the  second  of  equations  (3.5)  implies  that 


(3.6) 


rs  - 


u , (s , T|)  = - (1  - e'b)v.  01)  + e | eV(5,TDdg  ; 
1 l i 


hence,  u^O,!])  = 0 and  u^s,!)) 


s=0 


= - v (TD  for  all  T|£n.  Defining 


(3.7) 


a s rm(s)(l-e  s)ds  , 

«t  _ , 

( Tl)  S w.(Tl)  - J ra(s)e 


e'z^  (g.Tl)dg  ds  , 


we  note  that  w£R(ri)  H (C  (fl))^,  and  (2.3)  implies  that  0 < a < ».  Using 
(3.6)  and  (3.7),  the  first  of  equations  (3.5)  becomes 


v.(T\)  - 


(3.8) 
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It  follows  from  elliptic  theory  (see  Mizohata  (1973))  that  (3.8)  admits  a 
solution  v€R(Il)n(C  (Q))^  such  that  v.(Tl)  = 0 for  11 6 T;  inserting  this  v 

oo  — 3 

in  (3.6)  we  see  that  u € (C  ([0,r)  x Q))  . If  r < ®,  we  now  can  conclude  that 
(v,u)  €£(A).  If  r = ®,  let  0 < 3 < ® and  note  that,  by  (3.5)  and  through 
the  same  argument  used  in  proving  that  -A  was  accretive. 


li(w»z)|lY  s {w  (T]) w (Tl)  + I m(s)3  z (s.Tl)  *3.z.  (s.TDds)  dQ 
x i i j i j i J 


[ iv.(7Dv.(T])  + ! m(s)3  .u.  (s,Tl) -3  .u.  (s.Tl)ds^  dO 

Jti  l Jq  ji  j i ^ 

+ ]^m(s)ajux(s^)ds'TTik3h  j^(s)Vk(s^ds 

+ j^m(s)3j  [v.(Tl)  + u.(s,Tl)]  *3.  [v.(TD  + ^ u.  (s , T»]  ds  } df 


Letting  3 - r = ®,  we  see  that  (u,v)  €f(A),  with  ||(u,v)||£  + ||a(u,v)||£  * ||(w,z)]|£ 
Hence,  for  any  r ^ (0,®],  we  have  shown  that  (3.5)  has  a solution  (v,u)  € £(A) 
for  every  (w,z)  €S,  and  thus  we  conclude  that  R(I  - A)  = X.  The  proof  is  complete 


Theorem  3.1  shows  that  our  abstract  linear  evolution  equation  (3.4)  has 
a unique  solution  for  every  6 fi(A) . Furthermore,  as  the  dynamical  system 
lT(t)  )t  2-  o is  linear,  this  theorem  implies  that  all  positive  orbits  are  bounded 
and  every  motion  is  stable. 

For  r < ® the  injection  (£(A) , || • ||^)C^X  is  compact,  and  it  follows  that 
(I -A)  is  a compact  operator.  Hence,  if  r < ®,  we  may  now  conclude  that 
all  positive  orbits  are  precompact  (see  Dafermos  and  Slemrod  (1973)),  and 
this  fact  enables  us  to  prove  the  following  result. 
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Theorem  3.2:  For  m(s)  satisfying  (2.3)  with  r < ®,  every  motion  of 
is  asymptotically  stable. 


Proof:  Defining  V : X -•  ft  as  V(x)  = (x,x),  and  defining  V : X - R by 


V(x)  = lim  inf  - [v(T(t)x)  -V(x)]  , x€X 

t>kO  1 


we  see  that  V(x)  = 2(x,Ax)  for  x€t(A);  therefore,  our  accretivity  argument  shows 
that 


r <*r 

V(x)  = I 1 m'(s)*du(s,T\)-du.(s,Tl)dsdft  3 - W(x)  3 0 

Vo  k 1 k 1 


for  all  x * (v,u)  €£(A).  By  Theorem  3.9  of  Walker  (1976),  it  follows  that 

V(x)  £-W(x)s  0 for  all  x£X.  Hence,  V is  a Liapunov  function  on  X (see  Hale  (1969)). 

Conditions  (2.3)  on  m(s)  imply  that  m' (s)  < 0 for  all  s in  some  nonempty 
open  set  «$C(o,r);  consequently, 


[x€X|v(x) 


(s,td  •akui 


(s , T|)  dfi  = 


0 


a.e . 


Using  equation  (3.4),  it  is  not  difficult  to  see  that  the  largest  positive 
invariant  subset  of  {x€X|v(x)  * 0}  is  Ut+  * [o}.  As  all  positive  orbits 
are  precompact,  LaSalle's  Invariance  Principle  (see  Hale  (1969))  now  implies 
that  T(t)x  -»  0 as  t -»  ®,  for  every  x€X;  hence,  x ■ 0 is  an  asymptotically 
stable  equilibrium.  By  the  linearity  of  the  dynamical  system,  it  follows  that 
all  motions  are  asymptotically  stable,  and  the  proof  is  complete. 


In  two  recent  papers  Slemrod  (1976),  (1978)  has  used  similar  methods  to 
study  a problem  very  closely  related  to  our  formal  linear  problem  (3.1),  (3.2), 
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(2.14).  Rather  than  using  the  pair  (v,u)  to  represent  the  state  of  the  system, 

t a \ 

Slemrod  chose  the  pair  \v , - uj.  He  also  placed  additional  restrictions 
on  the  history  corresponding  to  (2.16)  and  a linearization  of  (2.17),  namely 

u(s,TD  + T|€ft  t (s.Tj)  €R+  x n , 

(3.9)  u(s,TD  - 0 , (s.Tj)  €R+  x p , 

aiui(s,Ti)  - 0 , (s,Tl)€R+xn  . 

These  conditions  require  the  "linearized  fluid"  to  have  a lways  occupied  the 
container  and  to  have  always  been  incompressible.  Using  a different  space 
and  topology,  Slemrod  (1976)  was  then  able  to  prove  that  his  formal  equations 
led  to  a stable  linear  dynamical  system,  paralleling  our  Theorem  3.1.  He 
also  obtained  an  asymptotic  stability  result  under  the  additional  assumption 

rr  2 

that  s m(s)ds  < ®,  which  does  not  presume  a finite  memory  length  r.  Our 
J0 

Theorem  3.2  assumes  r < <°,  but  does  not  require  the  additional  assumption  (3.9). 

Under  further  assumptions  on  the  behavior  of  the  material  function  m(s) , 
Slemrod  (1978)  proved  an  exponential  stability  result.  We  will  now  present 
a result  on  exponential  decay,  in  our  topology,  for  those  initial 
states  in  X that  happen  to  satisfy  (3.9).  To  this  end,  let  Q denote  the 
11*11^  - completion  of  the  set  t (v,u)  6 S|  (3.9)  holds}. 

Theorem  3.3:  Let  m(s)  satisfy  (2.3)  and  let  there  exist  ^ 2 ^ > ^ such  that 

5jm(s)  £ - m'(s)  * &2m(s)  *or  s € [o,r) 

Then  there  exist  M > 0,  e > 0,  such  that  ||T(t)x||j  * Me"et||x||1  for  all 
t € R+ , x € Q . 


IT 


17 


Proof:  It  is  not  difficult  to  see  that  4 is  positive  invariant  under 

t,T ( t ) } 2 g ; i.e.,  Y(x)  C-C}  for  x 6 (i.  Defining  U : X -•  R as 

U(x)  s j|xi|3  + 0 I v.  (T|)  | 1:1(3)0.(5,11)03  00  , x€X  , 

X 1 o0  i 


where  0 > 0,  we  fin0  that,  for  x€Q, 


U(x)  s lim  inf  -p  [U(T(t)x)  -U(x)] 
t\»0 


= - aS  ' v.  (T\)v.  (Tl)dO  + B [ v.(T0  \ m'  (s)u.  (s,Tl)0s  00 
J0  1 1 c0  1 J0  1 

+ 0 !"  ( 1 m(s)S  .u.  (s,Tl)dsV  , 10(5)0.0.(5,11)03)  00 

V "0  J 1 ' ' J0  J 1 

, fr ' 

+ j j m' (s)3^oi(s,TD  •0yai(s,Tl)0s  00  , 


where  a = 1 m(s)0s. 

J0 

We  wish  to  show  the  existence  of  numbers  c^  2 c^  > 0,  e > 0,  such  that 
cl!|x||£  2 u (x)  ^ c2||x||2  an0  U(x)  £ - 2eU(x)  for  all  x€Q.  We  first  notice 
that,  by  Schwarz'  inequality, 


m(s)0  u (s,Tl)ds  * a , m(s)3  u (s,Tl)  *0  u (s,T\)ds  . 

0 J 1 "0  ^ 1 J 1 


Also,  for  all  z € (C  (0))  it  is  known  that 
o 


[ z (T»z  (Tl)dO  s k [ 0 z (Tl)-0  z (TD0O 
1 1 Jn  J 1 J 1 


for  some  k(0)  < ®;  hence,  for  (v,u)  €Q, 


r i* 

| v.  (7))  I n'  (s)u,  (s  , TJ)  ds  <tt 
1 v 0 i 


s I f f m'  (s)v.  (71)  v.  (TD ds  dfi 

•Vo  1 1 


!*  f 

I m' (s)u . (s , T\)u  (s , T\) ds  dfi 

1 a o 1 1 


m(0)k  !*  v (71) v (71) dfl*  f f m' (s)d,u  (s,Tl)  u (s,71)ds  dfl  | 

1 1 JnJo  j 1 J 1 ' 

m(0)k?9  f v.  (71)v.  (TDdfl*  f f m(s)  3 u (s  , T\)  • 3 .u  (s , 71)  ds  dfl  , 

2 Jfl  1 1 "ajo  J 1 j i 


where  we  have  applied  Schwarz'  inequality.  Similarly)  for  x€Q, 

2 


» «r 

v.  (Tl)  m(s)u  (s.TQds  dfl 
' Q 1 JQ  X 


i ak  1 v.  (Tl)v  (ID dfl.  S'  f m(s)3  u ( s , T\)  -3  u (s,T])ds  dfi 
Jfl  1 1 ‘ CT  0 j 1 j 1 


It  then  follows  that 

]u(x)  - ||xi|^|  s (3ak/2)||x||£  , 


U(x)  s - aS  ; v (T\) v (71)  dfl 
VQ  1 1 


B(kl2m(0))^'^  j*  vi(TDvt(Tl) dfl)  ^ mCs^i^  (s  ,71)  "d^u^s  ,71)ds  dfl  J 


rr 


- (?,  -«P)'  m(s)d.u  ( S , H) -3  u (s.TDds  dfi 

1 Vo  j 1 j 1 


2 

for  all  x6^.  Choosing  3 > 0 so  small  that  3ak  < 2 and  3(k£2m(0)  + 4a  ) < 4a(^, 
we  see  that  suitable  numbers  c^,  c^,  € do  exist.  As  4 is  positive  invariant, 

.Off  J. 

it  follows  that  U(T(t)x)  * e U(x)  for  all  t£ft  , x£Q;  hence,  we  find  that 


||T(t)x||£  S (c1/c2)e'2£t|ix||^ 


t€R+  , x € Q , 


and  the  proof  is  complete. 
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In  this  section  our  principal  purpose  was  to  prove  that  the  linear  abstract  < 

evolution  equation  (3.4)  generates  a dynamical  system  on  X.  We  have  gone 
beyond  this  objective,  considering  stability  properties  and  exponential  decay 
of  motions,  in  order  to  provide  a basis  for  comparison  with  the  related 
results  of  Slemrod  (1976),  (1978). 


i 
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4.  The  Nonlinear  Problem 

We  now  return  to  the  nonlinear  problem  described  by  the  abstract  evolution 
equation  (2 . 19)  , 

(4.1)  x(t)  = Nx(t)  a.e.  t€ft+  , 

x(0)  = x € X , 

o 

where  N is  the  closed,  densely  defined  operator  described  in  §2.  Questions 
regarding  existence  and  uniqueness  of  solutions  of  (4.1)  are  quite  difficult 
to  resolve,  and  we  are  not  able  to  prove  that  (4.1)  is  associated  with  a 
dynamical  system  on  X.  However,  assuming  that  it  is,  we  shall  show  in  this 
section  that  all  positive  orbits  are  bounded  and  the  equilibrium  at  x = 0 
is  stable. 

In  §3  we  showed  that  linearization  of  the  problem  led  to  a dynamical 
system  on  the  Hilbert  space  X.  This  suggests  that  it  is  plausible  to  assume 
that  the  nonlinear  problem  also  is  associated  with  a dynamical  system  on  X, 
in  the  following  sense. 

Assumption  4.1:  For  all  sufficiently  small  \ > 0,  ft  (I  - MO  = X and  = (I  - MO 
exists;  moreover,  with  s I,  there  exists  a dynamical  system  [s(t)}t  on  ^ 
such  that  j”,  x - S(t)x  as  n -*  ®,  for  every  x 6 X,  t €ft+,  the  convergence  being 
uniform  on  compact  subsets  of  R+. 

This  particular  association  between  the  dynamical  system  {.S(t)  ^ ^ and 

the  operator  N is  motivated  by  considering  a backward-difference  approximation 
of  equation  (4.1)  given  by 
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Ant'  Ant  t\  t Ant  > , r B+ 

\n  1 \n  n/  n \jn  / ’ » » . » » 


x(0)  = x 

o 


It  is  seen  that  for  sufficiently  large  n,  depending  on  t,  this  equation  has  a 

solution  x(t)  = Jn.  x if  ft(I  - XN)  = X and  if  J.  = (I  - XN)  1 exists  for  all 
t/n  o \ 

sufficiently  small  X > 0.  Under  Assumption  4.1,  J^nx0  is  a Hille-type 
approximation  of  S(t)xo  (see  Yosida  (1978)). 

If  0)1  - N were  accretive  for  some  U)€ft,  then  the  theory  of  Crandall  and 
Liggett  (1971)  would  show  that  Assumption  4.1  holds  if  and  only  if  ft(I  - XN)  = X 
for  ail  sufficiently  small  X>  0 (see  Yosida  (1978)).  Unfortunately,  u)I  - N 
is  not  accretive  and  we  are  unable  to  prove  the  validity  of  our  assumption, 
even  if  ft (I  - XN)  = X, 

It  is  remarkable  that,  under  Assumption  4.1,  it  is  easy  to  prove  that 
all  positive  orbits  are  bounded  and  the  equilibrium  at  x = 0 is  stable. 


Theorem  4,2:  For  m(s)  satisfying  conditions  (2.3),  and  under  Assumption  4.1, 
|js(t)x||^  £ j|x|l^  ^or  t^^+>  x € X.  Furthermore,  S(t)x-  0 as  t - ® if 
the  positive  orbit  Y(x)  is  precompact. 

Proof.  Defining  V(x)  s ||x||^  and 

V(x)  s lin,  inf  i [v(S(t)x)  -V(x)J 

t VkO 

for  x€X,  we  note  that  for  X > 0 and  x€£(N), 
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V(x  - \Nx)  * V(x)  - 2\<x,Nx>  , • 

r r 'r  . -I 

(x,Nx>  = - j v.  ("H)  *'niX\Vj  (^)  *ajVX^>  + J m(s)  ^xUk^S’^^  ’SjUk<'S’^  + ajUX^s»^  Jdsj 

- f | m(s)d.u.(s,TD-a .{v.(Tl)  + v (Tl).S  u (s.TD  + §-  u (s.Ttf}  ds  dO 

•‘ffo  k 1 ^ J J 

* - J vj  (°)  * dj{v^ (T^)vj& (T\)  + j m(s)Skui(s,1ll)  •&kui(s,T\)dsj  dfi 

J* 

- [ v,(Tl)-3.  | m(s)[a(,uv(s,TD-a.uv(s>T\)  + , T\)  ]ds  dfi 


V*  j V * k J k 


. r 


[ a V (T\)-[  m(s)[a.u  (s.Tl) -S  u (s.Tl)  + 3 u (s , T|) ]ds  dfi 
K J -Q  J J 


• i InJom(s)  ti  c\vs  • ^ ’ Vi (s>  71)3(18  d0 

J m'  (s)  Ws,71)  •3ku^(s,T\)ds  dfi  = - W(x)  £ 0 


1 

2 J 


As  i>(N)  is  dense  and  Assumption  4.1  is  made,  Theorem  3.4  of  Walker  (1979) 
shows  that  V(x)  £ - W(x)  s 0 for  all  x € X;  consequently,  V(S(t)x)  s V(x)  for 
all  t£R+,  x € X.  As  in  the  proof  of  Theorem  3.2,  we  also  see  that  the  largest 
positive  invariant  set  m+  in  {x£X]v(x)  = 0}  is  rn+  = £o } ; hence,  if  x is 
such  that  Y(x)  is  precompact,  LaSalle's  Invariance  Principle  (see  Hale  (1969)) 
implies  that  S(t)x  -*  0 as  t -*  ®.  The  proof  is  complete. 


Theorem  4.2  shows  that  all  positive  orbits  are  bounded  and  the  equilibrium 
at  x * 0 is  stable.  If  all  positive  orbits  could  be  shown  to  be  precompact, 
Theorem  4.2  would  imply  that  x = 0 is  globally  asymptotically  stable. 

Although  we  do  not  know  that  all  positive  orbits  are  precompact,  the  last  result 
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of  Theorem  4.2  does  show  that  there  exist  no  equilibria  (steady  flows)  other 
than  x = 0 (the  rest  condition),  and  it  also  shows  that  there  are  no  nontrivial 
periodic  motions  (nonsteady  periodic  flows). 

VJe  remark  that  the  function  V (x)  s ||x|i^  used  in  the  proof  of  the  foregoing 
theorem  is  a Liapunov  function  for  [s(t)}t2g  (see  Ha^-e  (1969)  and  Walker  (1979) 
Useful  Liapunov  functions  often  are  extremely  difficult  to  discover  for  a 
highly  nonlinear  problem,  and  discovery  of  a topology  suitable  for  a state 
space  may  be  even  more  difficult.  In  fact,  these  difficulties  are  so  inter- 
related (see  Walker  (1976))  that  a "formal"  Liapunov  function  (for  the  formal 
equation)  often  is  sought  a priori,  as  a means  of  suggesting  a suitable 
topology  for  the  state  space.  This  is  what  led  us  to  set  equation  (4.1)  ar.d 
Assumption  4.1  in  the  particular  Hilbert  space  X,  rather  than  in  any  other 
metric  space. 

The  function  V is  the  only  Liapunov  function  that  we  have  been  able  to 
find  for  the  nonlinear  problem.  This  is  in  contrast  with  the  linearized 
problem  of  §2,  which  admits  an  infinite  family  of  useful  Liapunov  functions, 
and  leads  to  a corresponding  family  of  linear  dynamical  systems  on  state 
spaces  differing  in  their  topologies.  The  linear  dynamical  system  and 
Liapunov  function  of  Slemrod  (1976)  belong  to  this  family,  which  can 
be  defined  in  cerms  of  the  set  of  linear  operators  that  commute  with  the 
linear  operator  A of  (3.3).  However,  among  all  of  the  functions  in  this 
family,  it  appears  that  only  the  function  V employed  here  is  useful  with  the 
original  nonlinear  problem. 

Our  Liapunov  function  V admits  a simple  physical  interpretation.  To 
demonstrate  this  point  most  clearly,  we  return  to  the  formal  problem  of  §2  and 
relax  certain  of  our  assumptions.  Rather  than  assuming  that  v (T],  t)  = 0 for 
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(T), t)  er  x R+,  let  us  assume  only  that  fluid  can  not  cross  the  boundary  T of  Cl; 
that  is,  v.  (Tl,t)n.  (Tl)  = 0 for  (T|,t)  €F  X ft+,  where  n(Tl)  denotes  the  unit 
outward  normal  to  P.  Let  us  also  consider  a general  body  force  field  f^CHjt), 
rather  than  one  derived  from  a potential  q(T],t).  If  we  retain  all  other 
assumptions  of  §2,  formal  computations  lead  to 

d r rr 

— V(x(t)  = j j m' (s)3kui(s,Tl,  t>dku.  (s,Tl,  t)ds  dF  + 2P(t) 

i 2P(t)  , t6R+  , 

where  P(t)  is  the  external  power, 

P(t)  s [ f . (T),  t)v . (T|,  t)  dfJ  + j n.  (Tl)s  . . (T],  t)v.  (Tl,  t)  dF 

"Q  1 1 1 1J  1 

Consequently,  we  see  that  V has  the  basic  property  of  a "mechanical  energy 

function"  for  the  fluid.  When  all  assumptions  of  §2  are  applied, 

P(t)  = 0 and  ^ V(x(t))  £ 0 for  all  t£R+. 

The  assumption  of  "fading  memory"  played  a crucial  role  in  our  analysis; 

however,  we  remark  that  conditions  (2.3)  can  be  slightly  weakened.  It  is  not 

difficult  to  see  that  we  do  not  actually  need  m€C^(ft+)  and,  except  for  Theorem 

all  of  our  theorems  continue  to  hold  if  conditions  (2.3)  are  replaced  by 

+ (* 

a)  m : R -*  fl  is  integrable,  with  0 < I m(s)ds  < “ , 

J0 

(4.2)  b)  m is  nonnegative  and  nonincreasing  on  R+  , 

c)  m(s)  > 0 for  s€[0,r)  ; m(s)  = 0 for  s£(0,r)  , 


where  r > 0 may  be  infinite. 
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